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We review the implementation of a model for multiple partonic interac- 
tions in HERWIG-I--I-. Moreover, we show how recent studies on the colour 
structure of events in Herwig-|— I- led to a significant improvement in the 
description of soft inclusive observables in pp interactions at the LHC. 
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1. Introduction 

Multiple partonic interaction (MPI) models are crucial for a successful 
description of the underlying event (UE) in hard hadronic collisions and 
minimum-bias (MB) data from the Tevatron and the Large Hadron Col- 
lider (LHC). The modern Monte Carlo event generators Herwig-|— |- [1], 
Pythia [2 [3] and Sherpa [4j all include an MPI model in order to simu- 
late the underlying event. In this short note we summarize the modelling 
of multiple parton interactions in the Herwig-|— |- event generator, present 
its recent developments and show a comparison of the improved model to 
LHC data. 

2. MPI model in HerwigH — |- 

The modelling of underlying events in Herwig-|— |- is based on the fact 
that at high enough energies the hard inclusive cross section for dijet pro- 
duction, 

^inc(^.pmin) = ^ / ^ dp^ f■/,^{x^, f,') ^ » fj/hA^2,^') , (1) 
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will eventually exceed the total cross section [5j. This leads to the in- 
terpretation that in fact the inclusive cross section counts not only single 
hard events but all hard events that occur in parallel during the very same 
hadron-hadron collision. 

2.1. Eikonal model 

With the eikonal assumption that at fixed impact parameter, b, indi- 
vidual scatterings are independenlj^ and that the distribution of partons in 
hadrons factorizes with respect to the b and x dependence, we can express 
the average number of hard interactions in a hadron-hadron collision as: 

where the function A{b; /i^) describes the spatial overlap of the two collid- 
ing hadrons (protons) as a function of the impact parameter b. In both 
Fortran Herwig [6J with a Jimmy plug-in [7j and Herwig-|-+, the over- 
lap function is modelled according to the electromagnetic form factor [8] 
and the parameter fi interpreted as the inverse radius of the proton. Since 
the spatial parton distribution (the colour distribution) is assumed to be 
similar to the distribution of electric charge, but not necessarily identical, 
fj, is treated as a free parameter. This model allows for the simulation of 
multiple interactions with perturbative scatters which have pt > p™™. 

The extension to soft scatterings (with pt < p™™) is kept as simple as 
possible. The additional soft contribution to the inclusive cross section is 
also eikonalized, such that we can also calculate an average number of soft 
scatters from the resulting fx^j^ and an overlap function Asoft{b) for the 

soft scattering centres. The functional form 74soft(&) is assumed to be the 
same as for the hard scatters, but we allow for a different inverse radius, 
/Xsoft- The consistency of this model with unitarity is given by fixing the 
two additional parameters cr™^^ and fJ-l^f^ from two constraints. First, we 
can calculate the total cross section from the eikonal model and fix it to be 
consistent with the Donnachie-Landshoff parametrization [9]. In addition, 
using the optical theorem, we can calculate the t-slope parameter from the 
eikonal model and fix it to a reasonable parametrization. This model and 
its extension to soft scatterings were implemented in Herwig-|--|- O |T0] 
and proved to be capable of describing the whole spectrum of UE data from 
the Tevatron [IH I12| (including its minimum bias part). However, the MPI 
model is too simple to describe the first MB data from ATLAS [1'6\ |T3] (see 
red line in Fig. [l]). A tuning of the parameters of the model |15j did not 



^ Of course some basic correlations are incorporated into the model later on, for exam- 
ple by requirement of momentum conservation. 
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Charged particle ij at yooGeV, track p_i_ > 500 MeV, for N^u > 6 Charged {p i_) vs. N^h at yooGcV, track pj^ > 500 McV, for Nd-, > 1 




Fig. 1. Comparison of Herwig 2.4.2 without CR and Herwig 2.5 with PGR to 
ATLAS minimum-bias distributions at ^/s = O.QTeV. 

improve the description. This observation led to new developments of the 
MPI model to include non-perturbative colour reconnections (CR) |16t I17j. 

2.2. Colour reconnection 

The colour structure of multiple interactions can cause non-trivial chang- 
es to the colour topology of the colliding system as a whole, with poten- 
tially major consequences for the particle multiplicity in the final state. The 
colour connections between partons define colour singlet objects known as 
clusters (see Fig. [2]). The cluster hadronization model [18], which is used 




Fig. 2. Formation of clusters, which we represent by ovals here. Colour lines are 
dashed. The right diagram shows a possible colour-reconnected state: the partons 
of the clusters A and B are arranged in new clusters, C and D. 

in Herwig-|— |-, generates clusters which are used as the starting point for 
the generation of hadrons via cluster decays. The idea of CR is based on 
colour preconfinement |19j, which implies that parton jets emerging from 
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different partonic interactions are colour-connected (clustered) if they are 
located closely in phase space. As the MPI model does not take that into 
account, those colour connections have to be adapted afterwards by means 
of a CR procedure pictured in Fig. [2] In our CR model we define the dis- 
tance between two partons to be small when their invariant mass (cluster 
mass) is small. Therefore, the aim of the CR model is to reduce the colour 
length A = J2i=i ) where Nd is the number of clusters in an event and mj 
is the invariant mass of cluster i. Based on this physical motivation, there 
are two colour reconnection models implemented in Herwig+-|-. The first 
CR model, so-called plain CR (PCR), has been included in Herwig++ as 
of version 2.5 [20]. This model iterates over all cluster pairs in a random 
order. Whenever a swap of colours is preferable, i.e. when the new cluster 
masses are smaller, this is done with a given probability (the only param- 
eter of the model). The PCR model has been shown to give the desired 
results (see blue line in Fig. [l]). However, as the clusters are presented to 
the model only in a given sequence, it is hard to assess which clusters are 
afi^ected and to what extent the sequence is physically relevant. Therefore, 
we implemented another model, the statistical CR model (SCR) [T7], which 
adopts the Metropolis algorithm to reduce A. This allowed us to systemati- 
cally study the consequences of the introduction of CR, the results of which 
are presented in details in [171 EI] ■ 



3. Tuning and Results 

Initially we were primarily aiming at an improved description of MB 
data, therefore we started by tuning the PCR model to ATLAS MB data[^ 
(results are shown in Fig. [T]). The next important question was whether 
the new model is able to describe the UE data sets at different collider 
energies. A dedicated tuning procedure [22] resulted in energy-extrapolated 
tunes UE-EE, for both PCR and SCR, in which all parameters are fixed 
except for p™"^, which varies with energy. These tunes allow us to describe 
data at different energies |1H (see Fig. [s]) with a simple parametrization 
of the p™™ energy dependence. This was an important step towards the 
understanding of the energy dependence of the model. Finally, we show the 
results of the model for energy flow as a function of rj for minimum-bias 
events at ^/s = 0.9 and 7 TeV (Fig. [4]) presented by the CMS collaboration 
during this workshop and also published in [25]. This observable was not 
available during the preparation of the tunes and these very good results 
can therefore be treated as a prediction of the model. The new model is 

Since currently there is no model for soft diffractive physics in Herwig++, we use the 
diffraction-reduced ATLAS MB measurement with an additional cut on the number 
of charged particles, Nch > 6. 
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Fig. 3. ATLAS data at 900GeV (1st column), CDF data at ISOOGeV (2nd column) 
and ATLAS data at 7TeV (3rd column), showing the multiplicity density of the 
charged particles in the "transverse" area as a function of p^'^'^. 
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Fig. 4. Energy flow as a function of rj for MB events at = 0.9 and 7 TeV. 



implemented and available with the recent tunes in Herwig++ version 2.6 



4. Summary 



We have summarized the multiple partonic interaction models in Her- 
WIG++ and expanded on the motivation and modelling of colour reconnec- 
tion. Furthermore, we have shown that (sufficiently diffraction-suppressed) 
MB data from the LHC and underlying-event observables are well described 
by the present model which makes it useful for the LHC collaborations. 
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